Derivation of Radiation from a Dipole immersed in a 
Homogeneous Magnetized Plasma 
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Abstract 

We are considering a dipole antenna immersed in cold plasma and investigate its radiation patterns. 
We are solving Maxwell's equations using (1) analytic solutions for the case when the antenna is parallel 
to the magnetic field Bo and (2) numerical solutions for any orientation and angle 9 of the antenna. 

We investigate various phase velocities regimes as well as the effect of varying the frequency and 
the plasma parameters on the electromagnetic waves are presented and discussed. 



1 Introduction 

In the 60's there seemed to be considerable interest 
QWB69] . |WB72j . [WB692] . [SHM] 1 exploring VLF 
radiation originating from dipolcs in cold, anisotropic 
homoegenous magnetized plasma. An anisotropic 
medium is important in the fabrication of modern 
radio electronics and in astrophysics and plasma 
physics applications among others. Furthermore, 
considering radiation from sources in a cold, uni- 
form magnetized plasma can provide insight into the 
problem of focusing of radiation [SHM . 

An especially interesting source to consider is 
a dipole antenna, as in plasma, it can launch EM 
waves in addition to electrostatic waves. Moreover, 



one outstanding one being a diagnosis of the plasma 
without disturbing it [SHM . Conventional diag- 
nostic devices such as probes for measuring elec- 
trostatic and magnetic fields not only contaminate 
the plasma but are too large for the investigation of 
the microscopic structure of the plasma. Thus this 
analysis would have two important scopes: a charac- 
terization of plasma and an exploration of antenna 
design, as efficiency would be quantified by the EM 
wave scattering by the plasma wave under investi- 
gation. 

Due to the complexity of parameter tensors in 
the dielectric tensor, solving Maxwell equations for 
an anisotropic magnetized plasma and deducing the 
radiation patterns for an electric dipole require in- 
the far field behavior of dipole antenna can be changed vo lved analytical calculations. The papers cited above 



by adjusting both the operating frequency and the 
applied magnetic field. Since no physical antenna 
radiates uniformly in all directions and it radiates 
EM waves better in certain directions than others, 
there is a need to understand radiation patterns. 
They possess a practical worth in providing insight 
into the coupling between VLF antenna systems and 
the magnetospheric plasma, as their details are es- 
pecially important if the satellite transmitter is a 
low-power device, since the possibility of generating 
waves of useful amplitude may depend on maximiz- 
ing the flow of power at a given wave normal dis- 
tribution |WB72j . |WB692j . The consideration of 
VLF emissions obtained using a transmitter aboard 
a magnetospheric satellite could provide informa- 
tion about optimizing the production of longitudi- 
nally propagating waves and increasing the output 
power of the satellite transmitter, thus of significant 
importance in the design of a VLF satellite trans- 
mitting system |WB69j . 

This study finds a variety of other applications, 



employed a variety of techniques such as saddle point 
integrations, Green function, full wave coordinates. 

The present paper considers a modern twist to 
the same type of problems, focusing on the radia- 
tion patterns arising from each of the spherical co- 
ordinates of the electromagnetic fields. We present 
our problem in Section [2] We obtain the fundamen- 
tal solutions of a system of Maxwell's equations for 
a cold, anisotropic, homogeneous magnetoplasma 
in Section |3.1| In particular, we apply these to a 
small antenna dipole in Section |3.2| and we con- 
sider the radiation patterns in Section [4] We then 
present several electromagnetic radiation patterns 
plots. Next, in Section [5] we investigate these with 
respect to various phase velocities, dipole lengths, 
frequencies and far-field radii at which the patterns 
are considered. We find interesting correlations of 
the lobe predominance to the increasing phase ve- 
locities and radii, as well as rotational symmetries. 
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2 Formulation 

2.1 Dielectric tensor and dispersion 
relations 

We solve Maxwell's equations for E and H outside 
of a dipole antenna in a magnetized plasma. For 
these derivations, we assume that the plasma is ho- 
mogeneous and characterized by a dielectric tensor 
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The relationship between the components of K 
and oj is given by: 
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where = a, = b, fcn = — and the w rP , w 
are the usual plasma frequencies. 

Using these ratios, we can write the components 
in the following way: 
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We add the complex parts to avoid singularities 
in the real plane, which should result in smoother 
the radiation patterns, and we vary the ratios of the 
plasma frequencies to explore various phase veloc- 
ities. Thus, we shall consider the ratios | as the 
varying phase velocities which will allow us to ex- 
plore the characteristics of radiation patterns. The 
reason varying phase velocities gives insight is that 
they correspond to various wave regimes as consid- 
ered in |ABB| . 

We will have the following expression for the di- 
electric tensor: 
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3 Derivations 



3.1 Analytic solutions for antenna par- 
allel to Bq 

We start with Maxwell's equations in cylindrical ge- 
ometry, in Fourier domain: 
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Using Bg for boundary conditions we find 
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where Hq is the Hankel function of zeroth kind, 
and we can solve for a and /3 by setting boundary 
conditions on E z and Bg at the antenna. 

3.2 Application to a small dipole 

We then consider a small dipole, where we have 

/ = I smf3(\z\ — L) 
V = V cos(3(\z\-L) (7) 

and we apply the Fourier transform with respect to 
z to get I s and V s . 

Matching Hg at an antenna of radius ro allows 
us to get 

E z (r ) = ai^ 1 Wo) + Z^o' Wo) 
E z (r ) = a Pl H { lY ( Pl r ) + l3p 2 H { Q lY (p 2 r ) 

B z (r ) = a C ^4 1 \ Pl r ) + p C -^H^(p 2 r®) 
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Now we can set up the following system: 
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using the dielectric tensor adds another level of com- 
plexity to the Maxwell equations to be solved. 
In that case the equations to be solved are: 
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where a = Z ^K ± - k 2 z cos 2 9, b = % - K x , c = 
(10) - fc 2 cos#sin#, C d = £ A- X) e = %K ± - k 2 z , f = 



Finally the electric field E z (r) will be given by: ^ fc 2 cos6 » sin6 ) ; g = i*^ _ k 2 sin9 
, N We then have 
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and the magnetic field B z (r) will be given by: 
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Similarly we can calculate the other components 
of the fields, E r , Eg, B r , Bg, using the formulas for 
EL and BL: 
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3.3 Varying the angle 6 between the 
antenna and Bq 

In the case where we vary the angle 9 between the 
antenna and the direction of the magnetic field B , 
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4 Radiation patterns 

A radiation pattern shows the relative far field strength 
of the quantity considered versus angular direction 
at a fixed distance from the antenna. Since we are 
dealing with a dipole antenna, the azimuth plane 
pattern is formed by slicing through a 3D pattern 
in the horizontal plane, the x — y plane in this case, 
while the elevation/polar plane pattern is formed by 
slicing the 3D pattern through an orthogonal plane 
(either x — z or y—z). One quantity we are especially 
interested in is the Poynting vector, which is defined 
as S ave = \Re{E x H*) = \Re{E^H* - E e H*). A 
plot of the time- averaged, Poynting vector, |iS a „ e | for 
a fixed distance r is also called the power pattern of 
the antenna, and it is important as it describes the 
total power radiated by an antenna. 

Here we show characteristic radiation patterns 
for various phase velocities for the Poynting vec- 
tor and each component of the Electric and Mag- 
netic fields. 



Polar pattern - pv vs theta 
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Figure 2: Polar Patterns for E, rsph = 500, L 
200, oj = 10000 



Figure 1: Polar Patterns for S, rsph = 200 and 
rsph = 500, L = 100, w = 10000 
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Polar pattern - bsr vs theta 
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Figure 3: Polar Patterns for H, rsph = 500, L = 
200, lu = 10000 

5 Analysis 

We begin by remarking a few interesting general 
characteristics that will allow us to draw connec- 
tions among various types of patterns and param- 
eters. In the VLF range in the magnetosphere, for 
frequencies higher than the proton gyrofrequency, it 



is known that in general the only propagating elec- 
tromagnetic mode is the whistler mode, and conse- 
quently only this mode can contribute to the radi- 
ation power from a VLF antenna. A main, center 
lobe due to a principal saddle point implies that al- 
most all the radiation output is contained in that 
area, and this produces a focusing effect. Other 
lobes result from radiation focusing that is due to 
the presence of inflection points on the whistler- 
mode refractive index surface. These side lobes can 
give rise to nulls which would correspond to re- 
gions in which the antenna does not transmit sig- 
nal/radiation. Thus, exploring how these depend 
on the antenna length, frequency, the angle it makes 
with Bq, and other plasma characteristics gives great 
insight in antenna design. 

A big lobe in the radiation pattern signifies that 
the antenna radiates energy in one main direction. 
Focusing effects due to whistler inflection points may 
also occur for ray directions which are inclined to the 
static magnetic field line direction. This would im- 
ply that the radiation patterns from a given source 
may possess not only one major lobe along the mag- 
netic field line direction, but also an inclined one at 
a finite angle. 

As described in |WB69| . the drop-off in wave 
energy from the lobe maximum can be rapid and 
thus the longitudinal focusing of VLF waves can be 
used as a means of obtaining enhanced wave ampli- 
tudes from a low-power satellite transmitter. The 
fact that some sources may give rise to radiation 
patterns with not only a major lobe along the mag- 
netic field line direction, but also another which is 
inclined at some finite angle to this direction, im- 
plies that the focusing effects may also occur for ray 
directions which are inclined to the static magnetic 
field line direction. Whether or not the longitudinal 
lobe will predominate in an experimental situation 
will depend upon the scale of homogeneity of the 
magnetospheric plasma as well as upon the details 
of the radiating source. 

We will investigate the radiation patterns as the 
phase velocity vary. We attempt to arrange the 
radiation plots in such a way to emphasize com- 
mon characteristics, with varying phase velocities, 
lengths, radii. 

The plots we are going to refer to in our special- 
ized analyses appear below: 
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Figure 4: First row: radiation patterns of E r at a spherical radius of 500, a frequency u> = 10000, and to 
an antenna of varying lengths: 10,50, 100,200. Second row: radiation patterns of the spherical Hg at a 
spherical radius of 500 and a frequency u> = 10000, and to an antenna of varying lengths: 10, 50, 100, 200. 
The third row: radiation patterns of H r at a spherical radius of 500 and a frequency uj — 10000, and 
to an antenna of varying lengths: 10,50,100,200. The fourth row: radiation patterns of the E r of an 
antenna of length 100, and considered at various spherical radius distances: 200, 300, 500, 1000. The fifth 
row: radiation patterns of the H r of an antenna of length 100, and considered at various spherical radius 
distances: 200, 300, 500, 1000. 
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In the following sections, we will consider the 
influence of various factors (phase velocity, antenna 
length, spherical radius distance, frequency) on the 
radiation patterns. These considerations will vary 
in importance depending on the particular applica- 
tion envisioned, for example, for the magnetosphcric 
satellite transmitter described in the introduction. 
If wave-particle interaction experiments involving 
particles trapped on field lines that intersect the 
satellite are of interest, then one needs to maxi- 
mize the radiation along Bq. However, if an experi- 
ment involving wave-propagation is needed, such as 
in |SA68] , then it is important to maximize the radi- 
ation at high wave-normal angles, and, particularly, 
use a parallel electric dipole antenna. 

As we will see through these analyses, the far- 
field behavior of a dipole antenna can be changed 
by adjusting the plasma frequencies, its length, ap- 
plied magnetic field, spherical radius distance. The 
angle that the antenna makes with Bq also plays an 
important role, and is the subject of future studies. 

5.1 Dependence on the phase veloc- 
ity 

The phase velocity corresponds to the ratio ^ : the 
electron gyrofrequency versus the electron plasma 
frequency. In other words, this quantity measures 
the angular frequency of the circular motion of an 
electron in the plan perpendicular to the magnetic 
field versus the frequency with which electrons os- 
cillate. [ABB details the regimes corresponding to 
various phase velocities. 

As observed in the plot above, the radiation pat- 
tern lobes vary significantly with increasing phase 
velocities. The power decreases with increasing the 
collision frequency w ce . A possible physical expla- 
nation for this fact would be that for infrequent 
collisions, the particles spiral out to large gyroradii 
where they fall through a large potential difference, 
therefore absorb more energy from the wave. 

5.2 Dependence on antenna length 

In a non-plasma medium, the radiation patterns get 
more interesting as the electrical length of the dipole 
antennas increase, as more side lobes appear for L > 
A. 

The plots corresponding to a varying antenna 
length are presented in the first and second rows for 
the spherical radial coordinate of the electric coordi- 
nate of the electric field and the elevation coordinate 
of the magnetic field for the following parameters: 
spherical radius distance r_sph 500, uj 10000, and 
lengths varying between 10 and 200. 

From the plots, we notice that as L increases, 



Polar pattern - esr vs theta 
^ 90° 




Figure 5: Radiation patterns of E r at r sp h = 1000, 
L = 200 and L = 100, ui = 10000 

the side lobes of the spherical coordinates of the 
electromagnetic waves smoothen out and enlarge, 
while the main lobe expands. This would imply 
that radiation is focused, and that the nulls of the 
antenna will enlarge. 

On the other hand, the lobes of the Poynting 
vector will expand as the length of the antenna in- 
creases. This would imply that a longer antenna is 
capable of emitting more radiation across a wider 
range of directions, which would reduce the size of 
the nulls. 

5.3 Spherical radial magnetic field, Bsr 

As noticed in the radiation patterns in the third row 
in Figurc[5] the characteristics of the patterns are 
peculiar for Bsr, as we notice that as L increases, 
the lobes ply to the sides. 

5.4 Dependence on frequencies 

We notice the relative distribution of radiated en- 
ergy is substantially independent of the frequency 
uj. 

As representative radiation pattern plots corre- 
sponding to a varying frequency for a dipole antenna 
we presented below the dependence of the spherical 
radial coordinate of the magnetic field for the fol- 
lowing situations: r sp/ ^ = 1000, L = 100, uj = 5000 
and uj = 10000. 
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Polar pattern - psr vs theta 




Polar pattern - psr vs tneta 




antenna length, frequencies (u>, co ce ), the angle that 
the antenna makes with Bq, the radius at which 
we are obtaining the radiation patterns. We have 
found that the lobes characteristics and predomi- 
nance vary substantially with increasing phase ve- 
locities, and they either expand or shrink depend- 
ing on the type of electromagnetic radiation pattern 
considered. Moreover, we have found that the distri- 
bution of the lobes depends on the antenna length, 
as expected in vacuum. In conclusion, our results 
show that radiation pattern reshaping is possible 
by varying the constitutive parameters or operating 
frequencies for the cold magnetized plasma. 

7 Acknowledgments 

I thank Pat Colestock and Mark Galassi of Los Alamos 
National Laboratory for discussions on the problem 
and software help. 



5.5 Dependence on the radius at which 
the field intensities are considered 

The plots corresponding to a varying antenna length 
are presented in rows four and five of Figure [5] They 
correspond to the spherical radial coordinates of the 
electric and magnetic fields for an antenna of length 
L = 100, u> = 10000 and various distances between 
200 and 1000. 

As the spherical radius distance at which the 
patterns are considered increases, we notice the lobes 
shrink and ply to the sides. This would imply that, 
as expected, radiation varies inversely proportion- 
ally with the spherical radius distance, and that the 
space between the lobes will increase, which means 
that the antenna nulls will enlarge, causing a loss of 
radiation /signal. 

6 Conclusion 

In this paper we have analyzed the problem of elec- 
tromagnetic radiation from sources in a cold magnc- 
toplasma through the use of cylindrical coordinates 
in Fourier domain. In applying the resultant formu- 
lation, we have considered a small antenna dipole 
oriented at various angles with respect to the mag- 
netic field. Assuming plasma parameters appropri- 
ate to the VLF range in the magnetosphere, we have 
obtained expressions for the electromagnetic waves 
and used these to compute the radiation patterns. 
We have then explored these through varying phase 
velocities, and considered their dependence on the 
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